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Abstract 

This study presents a method to predict the growth fluctuation of 
firms interdependent in a network economy. The risk of downward growth 
fluctuation of firms is calculated from the statistics on Japanese industry. 



1 Introduction 

Does an abrupt ill health of one firm have a big impact on the health of others? 

Many firms ended in disastrous failure during the worldwide financial crisis in 
2008. Since then, risk managers, executives, and investors have been strongly in- 
terested in the transmission of distress and the knock-on defaults between firms 
which are interdependent in a network economy |May 2010] , In this study, a 
model for such firms is formulated with stochastic differential equations. Prob- 
ability parameters on trades between firms can be inferred statistically, and the 
time evolution of the net-worth of the firms can be predicted. The conditional 
value at risk of the downward growth fluctuation of firms is calculated from the 
statistics on Japanese industry in 2005. 

2 Stochastic model 

A stochastic model for firms is presented. Time dependent variables cii(t) for 
i = 1, 2, • • • , N is the net-worth of the i-th firm at time t. The interplay between 
the firms governs the time evolution of di(t). The income by the sales to others 
lets ai(t) increase. Its rate of change is given by 4>ijQ'j{t) where the probability 
parameters </>y are constant. The expenditure on the purchases from others 
and labor wages lets a,i(t) decrease. The rate is — \i<ii(t) where Ai is constant. 
This is a special case (linear production function) of the model for the financial 
accelerator in credit networks [Delli Gatti 2010] . It is also similar to the model 
for evolutionary autocatalytic sets |Mehrotra~ 2009 . 

The stochasticity of the time evolution ensues from an unpredictably irreg- 
ular pattern of trades between firms. The number of trades obeys a Poisson 
distribution if the probability of a trade per unit time is constant. The ampli- 
tude of fluctuation is nealy the square root of the average. The time evolution 
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of di{t) is given by a system of stochastic differential equations in eq.([T]). The 
functional form of the Gaussian white noises ^\t) and £j B] (i) is not known. 
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Eq.([T]) is equivalent to the Fokker-Planck equation in eq.©. It is a partial 
differential equation, which describes the time evolution of the joint probability 
density function P(a,t) of probability variables dj at t. 
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The drift and diffusion coefficients in eq.© are given by eq.Q and ((4]). 
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Predicting ai(t), given <fi and A, is a forward problem. Inferring the value 
of 4> arL d A from the observation on aj(i) statistically is an inverse problem 
|Maeno 2010J . These problems are mixed under practical conditions. The values 
of some parameters are known, and some data on Oi(£) are given. Eq.Q is 
converted to a system of ordinary differential equations in eq. ([5]) and ([6]) , which 
describe the time evolution of the 1st and 2nd order moments /zj 1 \t) and ^ijy'(i). 
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Generally, the time evolution of the rn-th order moments is given by a system 
of linear differential equations in eq.([7]). The elements of the N m x 1 vector 
M H(t) are ^(t), mIT^W- , P&Ljv(*)- The TV™ x ^iV m matrix A^, 
and iV™ 1 x N" 1 ^ 1 matrix £^ m ' are calculated from and _B l?fe . 



d/X '^ 1(t) = AM^W^+B^VK 1 ]^). (7) 

The solution of eq.© with the initial conditions ^ m l(0) is given by eq.©. 

H [m] {t) = exp(A^ m h){ f expi-A^h^B^fi^-^it^dt' + fi [m] {0))- (8) 
Jo 
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Figure 1: Risk R\,'. '(t) of the representative firms as a function of t. The j-th 
sector is Transportation equipment. The i-th sector is one of 12 selected sectors. 



Approximately, P(a, t) is a multi-variate normal distribution with the mean 
and covariance /i^(t). The exact formula for P(a,t) is obtained by the 
Edgeworth series. It is an asymptotic expansion of P(a, t) in terms of cumulants. 
The logarithmic likelihood function L(<p, A) is obtained immediately, once eq.© 
is solved, given a dataset a^(td) at td for the observations d = 1,2, • ■ • . It is 
given by eq.@. The estimators <j> and A are those which maximize L. 

L{4>,\) = Y J ^g{P(a^\t d )M<t>A))- (9) 
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3 Growth fluctuation 

The risk of firms are defined, and calculated with the Leontief coefficients 
in the input-output model for N = 34 Japanese industry sectors in 200£|. 
The q quantile value at risk vj q ^ (t) of the j-th firm at t is defined by eq. (fit)]) 
|Martine~J aramillo 2010] . It is the net- worth at which the cumulative density 
is q. Pyi{aj,t) is the marginal probability density function of dj. 

A j q] (*) 

P M (a j ,t)da j = q. (10) 

X) 

The conditional value at risk cj?j (t) which the ill health of the j-th firm 
imposes on the i-th firm at t is defined by eq. ([TTj) . Pc{a-i, t\ a j) IS the probability 



density function of conditioned on the value of cij ■ Note cj?j (t) ^ C|?l (t) . 



Pc(a i ,t|K [9l W)dai = g- (11) 



1 Ministry of Internal Affairs and Communications, 2005 Input-Output Tables for Japan, 
http: / /www. stat.go.jp/english/data/io/io05. htm 
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The quantity (t) is the risk of downward fluctuation from the expected 
growth at t. If P(a,t) is a multi-variate normal distribution, iiM(t) is given by 

cq.(fT2]l where ip^ = erf _1 (2g — 1). For example, f/jl 05 ! = —1.16 for 5 percentile 
and ^I 01 ! = —1.65 for 1 percentile. 



R w ( t ) = ihl ' ^ K ' = V — \/ ^— — v M ■ (12) 

The Leontief coefficients determine </>. The past growth rates of the sectors 
are used to obtain A. Suppose a representative firm in each industry sector 
whose share of production is 1%. Figure [T] shows i?^ 01 ' (£) of the firms in 
12 selected sectors (i) as a function of t. The j-th sector is Transportation 
equipment. The risk increases as time goes by. The representative firm in 
Mining have the largest risk, R « —0.6 at a quarter later (t = 90 days), when 
the representative firm in Transportation equipment falls ill. It is followed by 
the firms in Office supplies, Textile products, Non-ferrous metals, Electronic 
parts, Finance and insurance, Precision instruments, Iron and steel, Chemical 
products, Commerce, and Public administration. The firm in Medical service, 
health, social security and nursing care has the smallest risk. 
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